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ABSTRACT 

Current processing of acoustic array data is 
burdened with considerable uncertainty. This study 
reports an original methodology that serves to 
demystify array results, reduce misinterpretation, and 
accurately quantify position and strength of acoustic 
sources. Traditional array results represent noise 
sources that are convolved with array beamform 
response functions, which depend on array geometry, 
size (with respect to source position and distributions), 
and frequency. The Deconvolution Approach for the 
Mapping of Acoustic Sources (DAMAS) method 
removes beamforming characteristics from output 
presentations. A unique linear system of equations 
accounts for reciprocal influence at different locations 
over the array survey region. It makes no assumption 
beyond the traditional processing assumption of 
statistically independent noise sources. The full rank 
equations are solved with a new robust iterative 
method. DAMAS is quantitatively validated using 
archival data from a variety of prior high-lift airframe 
component noise studies, including flap edge/cove, 
trailing edge, leading edge, slat, and calibration 
sources. Presentations are explicit and straightforward, 
as the noise radiated from a region of interest is 
determined by simply summing the mean-squared 
values over that region. DAMAS can fully replace 
existing array processing and presentations 
methodology in most applications. It appears to 
dramatically increase the value of arrays to the field of 
experimental acoustics. 
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SYMBOLS 

a m shear layer refraction amplitude correction for 

L// 

A DAMAS matrix with A nn , components 
A nn > reciprocal influence of beamforming 

characteristics between grid points 
B array “beamwidth” of 3 dB down from beam 
peak maximum 

c 0 speed of sound without mean flow 

CSM cross spectral matrix 

D nominal diameter of array 

DR diagonal removal of G in array processing 
e steering vector for array to focus location 

e m component of e for microphone m 
f frequency 

A f frequency bandwidth resolution of spectra 

FFT Fast Fourier Transform 

< p array elevation angle 

G mm f cross-spectrum between p m and p m , 

G matrix (CSM) of cross-spectrum elements 

^ mm' 

H height of chosen scanning plane 

i iteration number 

k counting number of CSM averages, also 

acoustic wave number 

l representative dimension of source geometry 

detail 

LADA Large Aperture Directional Array 
LE leading edge 

m microphone identity number in array 
m' same as m , but independently varied 
m 0 total number of microphones in array 
n grid point number on scanning plane(s) 

M wind tunnel test Mach number 

N total number of grid points over scanning 
plane(s) 

p m pressure time records from microphone m 
P m Fourier Transform of p m 
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QFF Quiet Flow Facility 

Q„ idealized P m for modeled source at n for 
quiescent acoustic medium 
r c distance r m for m being the center c 
microphone 

r m retarded coordinate distance to m , r m c 0 

R nominal distance of array from scanning plane 

SADA Small Aperture Directional Array 

STD standard or classical array processing 

T complex transpose (superscript) 

TE trailing edge 

x m propagation time from grid point to 

microphone m 

w m frequency dependent shading (or weighting) 

for m 

W shading matrix of w m terms 

W width of scanning plane 

Ax widthwise spacing of grid points 

X matrix of X n terms 

X n “noise source’’ at grid point n with levels 
defined at array, Q*Q„ 

Ay heightwise spacing of grid points 

Y (e) output power response of the array at focus 

location 

Y matrix of Y n terms 

Y n Y ( e ) , when focused at grid point n 

Subscripts 

bkg background 

diag diagonal 

m : n term associated with m, as it relates to grid 
position Tl 
mod modeled 

INTRODUCTION 

Over recent years, particularly since the mid 
1990’s, the use of phased arrays of microphones in the 
study of aeroacoustic sources has greatly increased. Its 
popularity is due in large part to the seemingly magical 
presentations of array-processed results, which pull out 
features (real or just perceived) of noise source 
distributions on both wind tunnel models and full-scale 
aircraft. Properly utilized, arrays have been powerful 
tools that can often be used to extract noise source 
radiation information in circumstances where other 
measurement techniques may fail. However, 
presentations of array measurements of aeroacoustic 
noise sources can lend themselves to great uncertainty 
in interpretation. Proper interpretation requires 
knowledge of the principles of phased arrays and 
processing methodology. Even then, because of the 


complexity and differing processing philosophy, there 
can be misinterpretations of actual source distributions 
(and subsequent misdirection of engineering efforts). 

Before the mid 1980’s for aeroacoustic studies, 
processing of array microphone signals involved time 
delay shifting of signals and summing in order to 
strengthen contributions from, and thus “focus” on, 
chosen locations over surfaces or positions in the flow 
field. Over the years, with great advances in 
computers, this basic “delay and sum” processing 
approach has been replaced by “classical 
beamforming” approaches involving spectral 
processing to form cross spectral matrices (CSM) and 
phase shifting using increasingly large array element 
numbers. Such advances have greatly increased 
productivity and processing flexibility, but have not 
changed at all the interpretation complexity of the 
processed array results. Brooks and Humphreys 1 
examined current capabilities and determined the role 
of array size on the quantitative measurement of 
airframe noise sources. 

At Langley Research Center, for testing in the 
Quiet Flow Facility (QFF), quantitative definition of 
different airframe noise sources spectra and directivity 
has been a primary goal. This goal has been achieved 
with arrays in a rather straight forward manner for the 
localized intense source of flap edge noise 2 ; although 
for precise source localization. Coherent Output Power 
(COP) methods were used by incorporating unsteady 
surface pressure measurements along with the array. 
Quantitative measurements for distributed sources of 
slat noise 3 were achieved using an array and specially 
tailored weighting functions that matched array 
beampatterns with knowledge of the line source type 
distribution for slat noise. Similar measurements for 
distributed trailing edge noise 4,5 and leading edge noise 
(due in this case to grit boundary layer tripping) were 
made along with special COP methodologies involving 
microphone groups. 

Outside of Langley, there have been a number of 
efforts to analyze and develop more effective array 
processing methodologies in order to more readily 
extract source information. Several efforts include 
those to better account for array resolution, ray path 
coherence loss, and source distribution coherence 6 and 
for test rig reflections 7 . In a simulation study 8 of 
methods for improving array output, particularly for 
suppressing sidelobe contamination, three 
beamforming techniques were examined. These were a 
cross spectral matrix (CSM) element weighting 
approach, a robust adaptive beamforming, and the 
CLEAN algorithm. The CSM weighting approach 
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reduced sidelobes compared to classical beamforming 
with some overall improvement in main resolution. 
The results for the adaptive beamformer, used with a 
specific constant added to the CSM matrix diagonal to 
avoid instability problems, were encouraging. The 
CLEAN algorithm was found to have the best overall 
performance for the simulated beamforming exercise. 
CLEAN is a deconvolution technique first used in 
radio astronomy 9 . It was also examined in another 
study 10 , along with a related algorithm called 
RELAX 11 , using Langley QFF experimental array 
calibration data for a no-flow condition. It found mixed 
success in separating out sources. In a different study 12 , 
using the same data, two robust adaptive beamforming 
methods were examined and found to be capable of 
giving sharp beamwidths and low sidelobes. It should 
be mentioned that the above methods, although perhaps 
offering promise, have not been demonstrated to 
produce quantitatively accurate source amplitudes and 
distributions for real test cases. In the particular method 
of CLEAN, questions have been raised 13 with regard to 
practicality of the algorithm for arrays in reflective 
wind tunnel environments. A method that did show 
promise with wind tunnel aeroacoustic data is the 
Spectral Estimation Method 14 (SEM). The method 
requires that the measured CSM of the array be 
compared to a simulated CSM constructed by defining 
distributions of compact patches of sources (or source 
areas) over a chosen aeroacoustic region of interest. 
The difference between the two CSMs is minimized 
using a Conjugate Gradient Method. The application of 
positivity constraints on the source solutions was found 
to be difficult. The resultant source distributions for the 
airframe noise cases examined were regarded as being 
feasible and realistic, although not unique. 

The goal of the present effort was to develop a 
complete deconvolution approach for the mapping of 
acoustic sources (DAMAS 15 ) to demystify 2D and 3D 
array results, to reduce misinterpretation, and to more 
accurately quantify position and strength of 
aeroacoustic sources. Traditional presentations of array 
results involve mapping (contour plotting) of array 
output over spatial regions. These maps do not truly 
represent noise source distributions, but ones that are 
convolved with the array response functions, which 
depend on array geometry, size (with respect to source 
position and distributions), and frequency. The 
deconvolution methodology presented in this paper 
employs these processed results (array output at grid 
points) over the survey regions and the associated array 
beamforming characteristics (relating the reciprocal 
influence of the different grid point locations) over the 
same regions where the array’s outputs are measured. 
A linear system of “A” (number of grid points in 


region) equations and “A” unknowns is created. These 
equations are solved in a straight-forward iteration 
approach. The end result of this effort is a unique 
robust deconvolution algorithm designed to determine 
the “true” noise source distribution over an 
aeroacoustic source region to replace the classical 
beamformed distributions. Example applications 
include ideal point and line noise source cases, well as 
conformation with well documented experimental 
airframe noise studies of wing trailing and leading edge 
noise, slat noise, and flap edge / flap cove noise. 

DAMAS ANALYSIS 

Beamforming. The first step in the analysis is to 
beamform over the source region, using what have 
become traditional methods. Post processing of 
simultaneously acquired data from the microphones of 
an array begins with computation of the cross-spectral 
matrix for each test case data set. The computation of 
each element of the matrix is performed using Fast 
Fourier Transforms (FFT) of the original data 
ensemble. The transform pairs P m (f, T) andP m -(/,T) 
are formed from pressure time records p m (t) and 
p m < (t) , defined at discrete sampling times that are At 
apart, of data block lengths T from microphones m 
and m' , respectively. The cross-spectrum matrix 
element is 1 

K 

G m ,Af) =-^^ J [ p Uf^) p Af^)\ a) 

s± k=\ 

This one-sided cross-spectrum is averaged over K 
block averages. The total record length is T tot = AT . 
The term w s is a data-window (such as Hamming) 
weighting constant. G mm > (/) is seen to be a complex 
spectrum with values at discrete frequencies /, which 
are A f apart. The bandwidth is A/' = 1 / T (Hz). The 
full matrix is, with m 0 being the total number of 
microphones in the array, 



Note that the lower triangular elements are complex 
conjugates of the upper triangular elements. 

The cross-spectral matrix is employed in 
conventional beamforming approaches to electronically 
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FIGURE 1 Illustration of open jet test configuration where the 
array is out of flow and the scanning plane positioned over 
aeroacoustic source region. 

“steer” to chosen noise source locations about an 
aeroacoustic test model. Figure 1 illustrates a particular 
test setup of a distribution of microphones of a phased 
array located outside the flow field containing an 
aeroacoustic model. A scanning plane of grid points is 
defined over the noise source region. (A scanning plane 
may, for example, be placed through the chordline of 
an airfoil section when studying trailing edge and/or 
leading edge noise.) The beamforming approach 
involves steering vectors associated with each 
microphone with respect to the chosen steering 
location. In Fig. 1 , the steering location is designated as 
grid point n . The steering vector is 

e = col[ej e 2 - e m J (3) 

where the component for each microphone m is 

e m -a m — exp{./'2#T,„} (4) 

A 


contain terms that account for the mean amplitude and 
phase changes due to convected and refracted sound 
transmission through the shear layer to each 
microphone. The corrections are calculated 16 by the use 

17 

of Snell’s law in Amiet’s method , adapted to a curved 
three-dimensional mean shear layer surface defined in 
the shear layer. a m is the refraction amplitude 
correction. At mshear is the additional time (compared 
to a direct ray path with no flow) it takes an acoustic 
ray to travel to microphone m from the steering 
location n , due to the convection by the open jet flow 
and refraction by the shear layer. In Eq. (4), the ratio 
(r / r ) is included to normalize the distance related 

m C 

amplitude to that of the distance r. from the source 
location to the array center microphone at c. Both r 
and r c are in terms of “retarded” coordinates 516 . With 
this, r m =T m c 0 , where c 0 equals the speed of sound 
without mean flow. 

For classical or standard (STD) array 
beamforming, the output power spectrum (or response) 
of the array is obtained from 




e T Ge 

2~ 

m 0 


( 6 ) 


where the superscript T denotes a complex transpose 
of the steering vector. Flere Y[e ) is a mean-pressure- 
squared per frequency bandwidth quantity. The 
division by the number of array microphones squared 
serves to reference levels to that of an equivalent single 
microphone measurement. Note that the cross-spectral 
matrix (CSM) G often has a corresponding 
background cross-spectral matrix G bkg (obtained for a 
similar test condition except that the model is removed) 
subtracted from it to improve fidelity 1 . 

Shading algorithms can be used over distributions 
of array microphones to modify the output 
beampattern. The shaded steered response is 


The vector components serve to phase shift each 
microphone signal to allow constructive summing of 
contributions from the chosen locations. x m is the time 
required to propagate from grid point n to microphone 
m . In Ref. 1 . the phase is designated as 


T(e) 


e T WGW T e 

( '"o \ 2 

\ m =i ) 


( 7 ) 


m = (k-x m 'j + 2jfM 

m, shear ( 5 ) 

The term k is the acoustic wave vector, x m is the 
distance vector from the steering location to the 
microphone m. The steering vector components 


where w is the frequency dependent shading (or 

m 

weighting) for each microphone m . W is a row matrix 
containing the w terms. When all w terms are set to 

m m 

one and W becomes an identity matrix, all microphones 
are fully active in the beamforming to render Eq. (6). 
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(In Refs. 2-5 and 18,19, a special shading is used to 
maintain constant beamwidth over a range of 
frequencies by shading out ( w m = 0 ) inner microphone 
groups at low frequencies and by shading out outer 
groups at high frequencies.) 

A modified form of Eq. (6) is often used to 
improve dynamic range of the array results in poor 
signal-to-noise test applications. The primary intent is 
to remove the microphone self noise contamination 
(particularly caused by turbulence interacting with the 
microphones). This is done by removing (zeroing out) 
the diagonal terms of G and accounting for this change 
in the number of terms of G in the denominator. The 
output of Diagonal Removal (DR) processing is 


m 0 - - m 0 


DAMAS Inverse Problem Definition. The 

purpose here is to pose the array problem such that the 
desired quantities, the source strength distributions, are 
extracted cleanly from the beamforming array 
characteristics. First, the pressure transform P m of 
microphone m of Eq. (1) is related to a modeled source 
located at position n in the source field. 

Pnv.n (10) 

Here Q„ represents the pressure transform that P m . n 
(or P m ) would be if flow convection and shear layer 
refraction did not affect transmission of the noise to 
microphone m , and if m were at a distance of r c from 

n rather than r m . The e m . n ~ l term is simply those 
things that are postulated in Eq. (4) to affect the signal 
in the actual transmission to render P m . The product of 
pressure-transform terms of Eq. (1) becomes 


This modifies the beamform patterns compared to Eq. 
(6). The diagonal can be viewed as expendable 13 in the 
sense that it duplicates information contained in the 
cross terms of G . However, great care must be taken in 
physical interpretation of resulting array response maps 
- for example, negative “pressure-squared” values are 
to be expected over low-level noise source regions. The 
corresponding shaded version of Eq. (8) is 


p *p , = (O e~ l )* (O e~\ ) 

m:n m :n '^n^nv.n' 'X-nm :n > 


= n *0 (e~ l )*<?“* 


(ii) 


When this is substituted into Eq. (1), one obtains the 
modeled microphone array cross-spectral matrix for a 
single source located at n 


Y (.e) 


e T WG diag ^ 0 W T i 


( m 0 

1 

2 

1 m 0 ' 

\Zj W "' 

- 


\m=l 

) 


\m=l ) 


(9) 


^«mod X n 


( -1\* -1 
tei ) e x 

( -i\* -i 

te 2 ) 


/ -K* -i 

tei ) e 2 
( -i\* -i 

te 2 ) e 2 


/ -K * -1 

tel ) 


/ -1 x* -1 

(e m ) e m 
v m 0 m 0 


( 12 ) 


The common practice for studying aeroacoustic 
source of noise with arrays are to determine the array 
response, using either Eq. (6), (7), (8), or (9), over a 
range (grid) of steering locations about the source 
region. For particular frequencies, contours of the 
response levels are plotted over planes where sources 
are know to lie, or over volume regions in some cases. 
To extract quantitative contributions to the noise field 
from particular source locations, a number of methods 
are used. Integration methods, such as that of Ref. 1, 
can be used as well as special methods tailored to fit 
particular noise distributions, such as in Refs. 3-5. 
Still the methods can be difficult to apply and care 
must be taken in interpretation. This is because the 
processing of Eqs. (6)-(9) produces “source” maps 
which are as much a reflection of the array 
beamforming pattern characteristics as is the source 
distribution being measured. 


where X n is the mean square pressure per bandwidth at 
each microphone m normalized in level for a 
microphone at r m = r c . (A similar expression is used 
within the integration technique of Ref. 1.) It is now 
assumed that there are a number N of statistically 
independent sources, each at different n positions. One 
obtains for the total modeled cross-spectral matrix 



Employing this in Eq. (6), 




e T G 


mod c 


m n 


(13) 


(14) 
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where the bracketed term is that of Eq. (12). This can 
be shown to equal 




(16) 


where the components of matrix A are 


A' = - 




(17) 


By equating Y n (e) with processed F(e) from 
measured data, we have 


AX = Y 


( 18 ) 


Equation (18), for X, also applies for the cases of 
shaded standard, DR, and shaded DR beamforming, 
with components A nn > of A becoming 


e T W[ ] n ,W T e 

' m o \ 2 
\m = 1 / 


A' = - 


'([]„■) 


diag=0^n 


and 


A„„. = 


e T „W([ \ ,) diag ^W T e n 


' m 0 ^ 

2 

' >"o ' 


- 


\"'=1 ) 


\m = 1 / 


(19) 


( 20 ) 


( 21 ) 


respectively. For standard beamforming (shaded or not) 
the diagonal terms for A are equal to one. For 
Diagonal Removal beamforming (shaded or not), the 
diagonal terms for A are also equal to one, but the off- 
diagonal components differ and attain negative values 
when n and ri represent sufficiently distant points 
from one another, depending on frequency. 

Equation (18) represents a system of linear 
equations relating a spatial field of point locations, with 


beamformed array-output responses Y n , to equivalent 
source distributions X n at the same point locations. 
The same is true of Eq. (18) when Y n is the result of 
shaded and/or DR processing of the same acoustic 
field. X n is the same in both cases. (One is not 
restricted to these particular beamforming processing 
as long as A is appropriately defined.) Equation (18) 
with the appropriate A defines the DAMAS inverse 
problem. It is unique in that it or an equivalent equation 
must be the one utilized in order to disassociate the 
array itself from the sources being studied. Of course, 
the inverse problem must be solved in order to render 
X. 


DAMAS Inverse Problem Solution. Equation 
(18) is a system of linear equations. Matrix A is square 
(of size NxN) and if it were nonsingular (well- 
conditioned), the solution would simply be X = A~ X Y . 
However, it has been found for the present acoustic 
problems of interest that only for overly restricted 
resolution (distance between n grid points) or noise 
region size (spatial expanse of the N grid points) 
would A be nonsingular. Using Singular Value 
Decomposition (SVD) methodology 20 for determining 
the condition of A , it is found that for resolutions and 
region sizes of common interest in the noise source 
mapping problem in aeroacoustic testing that the rank 
of A can be quite low - often on the order of 0.25 and 
below. Rank here is defined as the number of linearly 
independent equations compared to the number of 
equations of Eq. (18), which is N= number of grid 
points. This means there are generally very large 
numbers of “solutions”. Equation (18) and the 
knowledge of the difficulty with equation rank were 
determined early in the present study. The SVD 
solution approach with and without regularization 
methodology 21 , special iterative solving methods such 
as Conjugate Gradient methods 20 , and others did not 
produce satisfactory results. Good results were 
ultimately obtained by a very simple tailored iterative 
method where a physically-necessary positivity 
constraint (making the problem deterministic) on the 
X components could be applied smoothly in the 
iteration. This is described below. 

A single linear equation component of Eq. (18) is 
A„ A + A n2 X 2 + . . . + A nn X n + . . . + A nN X N = Y n (22) 
With A nn = 1 , this is rearranged to give 
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n- 1 N 

X n =Y n -[^A nn ,X n , + 2)A m ,X n ,] (23) 

n'= 1 «'=«+l 

This equation is used in an iteration algorithm to obtain 
the source distribution X n for all n between 1 and N 
as per the following equation. 

N 

X ( l i) = Y 1 -[0 + 

n'= 1+1 

n- 1 N 

x i ° = - 1 2 a « X !? + 2 ] (24) 

n'=l n'=n+l 

jV-1 

^ ) =^v-[ 2 a ^' z «' )+o] 

n'=l 

For the first iteration (i = 1), the initial values X„can 
be taken as zero or Y n (the choice appears to cause 
little difference in convergence rates). It is seen that in 
the successive determination of X n , for increasing n, 
the values are continuously fed into the succeeding X n 
calculations. After each X n determination, if it is 
negative, its value is set to zero. Each iteration ( i ) is 
completed by like calculations, but reversed, moving 
from n = N back to n = 1 . The next iteration (i + 1) 
starts again at n = 1 . Equation (24) is the DAMAS 
inverse problem iterative solution. 

APPLICATION SIMULATION AND CRITERIA 

Key Application Parameters. Figure 2 has 
identified important parameters in defining the solution 
requirements for DAMAS for a scanning plane. The 
array has a spatial extent defined by the “diameter” D. 
It is at a nominal distance R from a scanning plane 
containing N grid points, which represent 
beamforming focal points, as well as the n locations of 
all the acoustic sources X n that influence the 
beamformed results Y n . For a particular frequency, the 
array’s beamformed output is shown projected on the 
plane as contour lines of constant output Y , in terms of 
dB. The scanning plane has a height of H and a width 
of W . The grid points are spaced Ax and Ay apart. 
Although not illustrated in Fig. 2, there are defined 
noise source sub-regions of size (. within the scanning 
plane (subsets of X n ), where details are desired. This 
relates to source resolution requirements and is 
considered below. For the scanning plane, the total 
number of grid points, 



FIGURE 2 Key geometric parameters of the array and source 
region scanning plane. 

N = [(W / Ax) + 1][(// / Ay) + 1] (25) 

The array beamwidth B is defined as the “diameter” of 
the 3 dB-down output of the array compared to that at 
the beamformed maximum response. For standard 
(STD) beamforming of Eq. (6), 

B ~ const x (R / fD) (26) 

For the SADA (Small Aperture Directional Array with 
a outer diameter of D = .65 feet) in a traditional QFF 
configuration 1 with R = 5 feet, the beamwidth 
isB *= (10 4 //) in feet for frequency / in Hertz. When 
using shading of Eq. (7), B is kept at about 1 ft. for 
10kHz</< 40kHz. 

In the applications of this report, some engineering 
choices are made with regard to what should represent 
meaningful solution requirements for DAMAS source 
definition calculations. Because the rank of matrix A 
of Eq. (18) equals one when using the iterative solution 
Eq. (24), there is no definitive limitation on the spacing 
or number of grid points or iterations to be used. The 
parameter ratios Ax IB (and Ay IB) and W IB (and 
H / B) appear to be most important for establishing 
resolution and spatial extent requirements of the 
scanning plane. The resolution Ax I B must be small or 
fine enough such that individual grid points along with 
other grid points represent a reasonable physical 
distribution of sources. However, too fine of a 
distribution would require substantial solution iterative 
times and then only give more detail than is 
realistically feasible, or believable, from a beampattern 
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which is too broad. On the other hand, too coarse of a 
distribution would render solutions of X which would 
reveal less detail than needed, and also which may be 
aliased (in analogy with FFT signal processing), with 
resulting false images. The spatial extent ratio W IB 
(and H / B) must be large enough to allow 
discrimination of mutual influence between the grid 
points. Because the total variation of level over the 
distance B is only 3 dB, it appears reasonable to 
require that 1 <W IB (and H/B). One could extend 
W IB (and H / B) substantially beyond one - such as 
to five or more. In the following simulations, resolution 
issues are examined for both a simple and a 
complicated noise source distribution. Two 
distributions types are considered because, as seen 
below with respect to £/B, source complexity affects 
source definition convergence. The simulations also 
serve as an introduction to the basic use of DAM AS. 

Regarding execution efficiency of the DAMAS 
technique, it is noted that the per-iteration execution 
time of the algorithm depends solely on the total 
number of grid points employed in the analysis and not 
on frequency-dependent parameters. In general, the 
iteration time can be expressed by time=C(2N) 2 i , 
where C is a hardware-dependent constant. A 
representative execution time is 0.38 seconds/iteration 
running a 2601 -point grid on a 2.8-GHz, Linux-based 
Pentium 4 machine using Intel Fortran to compile the 
code. For this study, a Beowulf cluster consisting of 
nine 2.8 GHz Pentium 4 machines was used to generate 
the figures shown subsequently. 

Simple Source. In a traditional contour type 
presentation, the top left frame of Fig. 3(a) shows an 
array output based on standard (STD) processing 
methodology of Eq. (6), being plotted in terms of 
constant dB contours over a scanning plane. In this 
simulation, the SADA is placed 5 feet from the plane 
that is positioned through a typical model location. In 
terms of the aforementioned parameters, H = W= 50” 
and Ax = Ay=l”. The resultant number of grid points is 
2601 (underlying grid points are not shown in top left 
frame). With a chosen frequency of 10 kHz and the 
beamforming of Eq. (6), B = 12”, so H / B=W / B= 
4.17 and Ax IB =0.083. A single synthetic point source 
is placed at a grid point in the center of the plane, at 
/7 = 1 30 1 . This is done by defining X 1301 to give 
100 dB = \QLogX nol and all other X values to zero in 

Eq. (18), and then solving for Y. The values of 
clB = 10 LogY n are then contour plotted. This, as with 
real array test data, is the starting point for the use of 
DAMAS. Equation (18) is solved for X using Eq. (17) 
for A nn ', by way of Eq. (24), using X n = Y n at the start 


Spanwise location (in) 

(a) /=10kHzand A.y/B=0.083 




DAMAS, 1000 Iterations 


-20 0 
Spanwise location (in) 

(b) /=20kHzand A.y/B=0.167 


"20 



Spanwise location (in) 
(c) /=30kHzand Ax!B= 0.25 


FIGURE 3 Synthetic point source - SADA STD processing for 
50”x50” scanning plane that is 5’ away. Resolution A.y=1”. 
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STD Beamforming 



Spanwise location (in) 
(a) /=10kHzand Aa7B=0.083 


STD Beamforming 



.DAMAS, 1000 Iterations 


-20 

Spanwise location (in) 
(b) /=20kHzand Ax/B=0.167 


20 


STD Beamforming 



DAMAS, 1 000 Iterations 


j".:. h jft. 


*!!* 


DAMAS, 5000 Iterations 


NASA 


0 20 -20 0 
Spanwise location (in) 

(c) /=30kHzand Ax/S=0.25 


20 


FIGURE 4 NASA image source - same conditions as Fig. 3. 


dB 60 70 80 90 100 



Spanwise location (in) 


FIGURE 5 Spatial aliasing with point source and image shifted 
between grid points. DAMAS with 1000 iterations. 

of the iteration. The bottom left frame of Fig. 3(a) 
shows the values of X n after one iteration (/' = 1) . 
Rather than showing contours, the presentation is one 
of X n values in terms of dB at the grid points. Each 
grid point is actually located at the bottom left corner 
of the “blocks”, each of dimension Ax = Ay = 1”. In 
the top right and bottom right frames of Fig. 3 (a), the 
results after the one thousandth (i = 1000) and the five 
thousandth (i = 5000) iteration, respectively, are 
shown. At the highest iteration value, the original input 
value of 100 dB has been recovered within 0.1 dB and 
that the surrounding grid values over the plane are 
down in level by about 40 dB, except for the adjoining 
grid points at about 15-20 dB down. At the lesser 
iteration numbers, although there is some spreading of 
the source region, the integrated (obtained by simple 
summing of values over the spread region) levels are 
very close to 100 dB. One obtains 99.06 dB for 100 
iteration (not shown in Fig. 3 (a)) and 100.03 dB for 
1000 iterations. 

The solution dependence on reducing the 
beamwidth B by a factor of two ( Ar/B=0.167) is 
demonstrated in Fig. 3(b) where the frequency used is 
20 kHz using the same standard processing over the 
same grid. The contour pattern is similar, but 
contracted, as shown in the left frame. The DAMAS 
result for 1000 iterations is given in the right frame. 
Comparing this to the results of Fig. 3(a), it is seen that 
here a more exact solution is attained with substantially 
less iterations. One obtains at peak of 99.97 dB, with 
all adjoining grid points lower by 27 dB. 
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The beamwidth 5 is again reduced for Fig. 3(c), 
where the frequency used is 30 kJfz, to obtain 
Ax IB =0.25. With the same result presentation as Fig. 
3(a), it is seen that again much less iteration is needed 
to attain more exact results with larger A xlB . At 1000 
iterations, one obtains 100.00 dB at the peak, with 
adjoining grid points being lower by 61 dB. Even at 
100 iterations (not shown), the results compare well 
with 5000 iterations of Fig. 3(a). 

Complicated Source. A more demanding 
simulation is shown in Fig. 4, where a distribution of n 
locations was defined with the sameX,, values (each 
corresponding to 100 dB) and others as zero. This 
gives a test of the solution procedure for a group of line 
source distributions. The scanning plane parameters, 
including the number of solution iterations, given for 
Fig. 4 are the same as those for Fig. 3. In Fig. 4(a) for 
Ax IB =0.083, the beamforming contour plot has an 
elongated appearance similar to one that would be 
expected for a single line source. But after using some 
iterations, one begins to see structure other than a 
single line source. Still, the image does not converge (it 
does not converge even after 30,000 iterations). In Fig. 
4(b), for Ax/5=0.167, a “very prominent” image 
emerges (the image is found to be recognizable even 
before 100 iterations). In Fig. 4(c), for Ax IB =0.25, all 
images are apparent. With regard to integrated power, 
it is found that when integrating (by summing grid 
point values), the total noise converges to the correct 
value rapidly with increasing iterations (within 1 dB 
after 16 iterations and within .05 dB after 100 
iterations) for all three Ax IB values. 

Comparing Figs. 3 and 4, it is seen that the degree 
of complexity of the source distributions does have an 
impact on the image convergence. It appears to be 
independent of the basic Ax/5 dependence 
demonstrated in Fig. 3. The source complexity 
resolution effect is quantified here by associating £, 
first, with a representative dimension of the closed 
portion of the letter “A” formed by the point sources (= 
6”) and, second, with the overall letter dimensions (= 
11”), and then relating these to beamwidth B. From 
this, and a subjective judgment of image quality, 
resolution of source region detail is only marginal for 
i / B < 1 , acceptable for 1 <,11 B , and good for 2< £ / B . 

Spatial Aliasing from Source Distribution. It is 

seen above that coarser Ax/5 resolutions require less 
iterations to get the same “accuracy”. And, this 
becomes even more true as the noise source becomes 
more complicated. A question is raised with regard to 
the potential error one may encounter if Ax IB is made 
too large in real data cases where significant sources 


may be in-between chosen grid points (all line sources 
contain energy between grid points). For any such 
error, analogy can be made with the common data 
analysis subject of aliasing errors with respect to FFT 
sampling rates. No problems of this nature are possible 
in the previous simulations because all sources are 
collocated at the grid points. To evaluate the degree of 
the potential aliasing problem, the simulations of Figs. 
3 and 4, are repeated with the noise sources being 
offset by one-half the distance between the grid points, 
thereby shifting the complete Y contour map 
diagonally by that amount. Figure 5 shows the 
DAMAS results for X n with this offset for two of the 
Ax/5 cases, one with Ax IB =0.25, corresponding to 
the coarsest cases of Figs. 3 and 4, and one twice as 
large at Ax IB =0.5. If there were no aliasing problem, 
the position offset for the single source should cause 
the four adjoining points to be equal and sum to a total 
of 100 dB, while the remaining scanning plane should 
have negligible summed levels. For Ax /5 =0.083, 
0.167, 0.25, 0.5, 0.75, and 1, the summed four levels 
for 1000 iterations (and the summed levels in the 
remaining plane) are 100.0 (67.8), 100.1 (82.7), 100.2 
(88.3), 100.5 (96.6), 99.8 (99.9), and 97.2 (97.5) dB, 
respectively. For the NASA image, as well as for the 
single source, it is found that there is small observable 
aliasing effects for Ax/5 =0.083 and 0.167. Flowever, 
it is seen in Fig. 5 that the NASA image brakes up to 
some degree for A x/5= 0.25, and to a larger extent for 
0.5. For Ax/5=0.75, the image is disrupted, and for 
Ax/5=1, it is destroyed. It should be kept in mind that 
these calculations represent limiting case tests for 
aliasing. Real data would be smoothed by the presence 
of other distributed sources. Still, it appears that to 
avoid aliasing problems, values of Ax/5 at or below 
0.2 are recommended. 

EXPERIMENTAL APPLICATIONS 

In this section, experimental data from several 
airframe component noise studies conducted in 
Langley’s QFF are re-examined with DAMAS. In 
these applications, DAMAS is not used with 
necessarily optimum resolution and scanning plane 
size. However, all cases fall at or near an acceptable 
range of 0.05 < Ax/5 (and Ay/5) < 0.2. For 
consistency with the simulations, (except for the 
calibrator case) the same scanning plane and resolution 
sizes are used with the same resultant number of grid 
points. The number of iterations used for all is 1000. In 
contrast with the simulations, the experimental results 
are presented in terms of one-third octave values, for 
the array using several different array beamforming 
methodologies, in order to compare to the results of the 
previous studies. 
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<|> = - 39 ° 


□ 
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(a) Noise path from the flap edge to the SADA is illustrated. 


Calibrator Source 



(b) Open end of calibrator source positioned next to flap 
edge. 

FIGURE 6 Sketches of test set up for flap edge noise test and 
calibration. 


shear layer refraction terms are absent in the steering 
vector definition, Eqs. (4) and (5). The right frame of 
Fig. 7(a) shows the result for the rendered source X 
distribution when DAMAS is applied, solving Eq. (18), 
using Eq. (19), by way of Eq. (24). The scanning plane 
used is H = W =12”, in order to match that of Fig. 7(a) 
of Ref. 1 . Consistent with the contour presentation, the 
DAMAS result is a one third octave presentation 
obtained by separately solving for the 546 separate 
bands and then summing. With B =12” and a chosen 
At = Ay =0.55”, one has a resolution of Ax / B =0.046 
(which is close to the recommended lower limit of 
0.05). The number of grid points is 441 and the number 
of iterations used is 1000 for each frequency. 


Note that a characteristic of the DAMAS solution 
is the non-negligible amplitudes distributed at grid 
points around the border of the scanning planes in Fig. 
7(b). This is a scanning plane “edge” effect that is 
found to occur only for experimental data, where noise 
in the scanning plane is influenced to some degree by 
sources outside (or extraneous to) the plane. DAMAS 
constructs noise distribution solutions on the scanning 
plane grid points totally based on whatever is measured 
by beamforming on those grid points. The edge effect 
was examined by expanding the scanning plane to 
eliminate any edge problem in the region of interest. A 
result almost identical to Fig. 7(a) was found over 
regions other than at the edge. Thus the edge effect has 
negligible impact on these results. This subject is dealt 
with subsequently for other applications. 


Calibrator Source w/wo Flow. A sketch of the 
flap edge noise experimental setup is shown in Fig. 
6(a). The airfoil main element is at 16° angle-of-attack 
to the vertical plane. The SADA is shown positioned 
out of the flow. For this configuration 1 , the calibration 
test is performed using a noise source, comprised of an 
open end of a one-inch diameter tube, placed next to 
the flap edge, as shown in Fig. 6(b). 

For the calibrator source operating with no tunnel 
flow M = 0, Fig. 7(a) shows SADA response contours 
for standard (STD) processing with shading, Eq. (7), 
over a scanning plane positioned through the airfoil 
chordline. This corresponds to a frame of Fig. 7 (a) in 
Ref. 1. This is a one-third octave presentation for 
10 LogY for / 1/3 = 40 kHz. The result was obtained by 
performing and summing 546 single-frequency 
beamforming maps (each with frequency resolution 
bandwidth of A f= 17.44 Hz). Note that with this array 
shading, only the inner SADA diameter of 1.95 inches 
is active. For this no-flow case, the convective and 


In Ref. 1, a small rectangular integration region, 
shown by dashed lines in Fig. 7(a), was used to 
calculate an integrated value of 62.8 dB. 
Correspondingly, for the present DAMAS result, one 
simply adds the pressure-squared values of the grid 
points within the source region. One obtains a value of 
62.9 dB. Figure 7(b) shows the SADA response 
contour for the tunnel flow at M= 0.17. Here, the 
convective and shear layer refraction terms are 
important in the steering vector definition. The 
integrated value from Ref. 1 is 58.1 dB, whereas the 
DAMAS value is 57.3 dB. It is seen by comparing the 
somewhat smeared image of Fig. 7(b) to 7(a) that the 
affect of the tunnel flow and the resultant turbulent 
shear layer is to spread the apparent noise region. The 
subject of turbulent shear layer noise scattering is 
analyzed in Refs. 1 and 22. The DAMAS result in Fig. 
7(b) is of particular interest because, to the knowledge 
of the authors, it may be the first direct measure of 
spatial dispersion of noise due to turbulence scatter. 
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(a) M=0, integrated level Ref.1 62.8 dB, DAMAS 62.9 dB 

STD Beamforming DAMAS, 1000 Iterations 



Spanwise location (in) 


(b) M=0.17, integrated level Ref.1 58.1 dB, DAMAS 57.3 dB 

FIGURE 7 Calibrator source test, SADA shaded STD 
processing, / 1/3 =40 kHz and Ax/.B=0.046 


For the same test cases as Fig. 7, Fig. 8 shows the 
results when diagonal removal (DR), Eq. (9), is 
employed in the beamforming. Correspondingly, 
DAMAS is applied using Eq. (21) for A nn f. It is seen 
that although the DR processing modifies the Y 
distributions, the X source distributions and values are 
calculated to be almost identical to those of Fig. 7. As 
mentioned previously, DR processing has the 
advantage of removing the autospectra (and possible 
microphone noise contamination) from the processing, 
while still maintaining full rank for the solution 
equations. 

Although it is beyond the scope of this paper to 
evaluate the use of DAMAS for different array designs 
than the SADA, a limited application using Large 
Aperture Directional Array (LADA) data produced 
good comparisons for a case corresponding to a frame 
of Fig. 7(a) in Ref. 1. The LADA has an outer diameter 
of D = 2.83 feet, which is 4.35 times the size of SADA 
(17.4 times the active diameter of the SADA at 40 kHz 
for shaded processing). For a similar calibration test to 
that of Fig. 7(a), but for array processing without 
shading, the integrated LADA value from Ref.1 is 60.3 
dB, and the corresponding DAMAS summed value is 
61.4 dB. 



DR Beamforming 


-5 0 5 

Spanwise 


DR DAMAS, 1000 Iterations 


-5 0 5 

location (in) 



(a) M=0, integrated level Ref.1 62.2 dB, DAMAS 61.6 dB 



(b) M=0.17, integrated level Ref.1 57.9 dB, DAMAS 56.9 dB 

FIGURE 8 Same conditions as Fig. 7, except shaded DR 
processing is used. 


Trailing Edge and Leading Edge Noise Test. 

The test configuration 4 where a NACA 63-216 airfoil, 
with a 16” chord and 36” span, is positioned at -1.2° 
angle-of-attack to the vertical flow is shown in Fig. 9. 
The flap is removed and the cove is filled-in such as to 
produce a spanwise uniform sharp Trailing Edge (TE) 
of 0.005”. Grit of size #90 is distributed over the first 
5% of the Leading Edge (LE) to ensure fully turbulent 
flow at the TE. The SADA position is at 0=90°. Figure 
10 shows the array output over a scanning plane placed 
through the chordline. The scanning plane of size 
H = W =50” extends “beyond” the sideplates that hold 
the 36” span airfoil. The sideplate regions as seen from 
the viewpoint of the array, to the left of the -18” 
spanwise location and to the right of the 18” location in 
Fig. 10, represent reflected source regions. 

The array output in Fig. 10 is presented for four 
one-third octave frequencies for STD processing with 
shading, Eq. (7). As before, individual frequency 
results are processed and are then summed to obtain the 
results shown. It is seen that for the / 1/3 = 3.15 kHz 
case, that the most intense region is just aft of the 
airfoil TE. As frequency is increased, the intense 
regions appear to first concentrate near the TE, and 
then to shift towards the LE. 
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FIGURE 9 Sketch of test set up for TE and LE noise test. 

Figure 1 1 presents DAMAS results corresponding 
to Fig. 10. For / 1/3 =3.15 kHz, Ax= Ay =1.8” is used to 
obtain the chosen lower-limit resolution of 
Ar/B=0.047. For 8, 12.5, and 20 kHz, the chosen 
Ar = Ay = l” give Ax / B=0.066, 0.083, and 0.083, 
respectively, for shaded beamforming. The results 
shown appear to very successfully reveal noise source 
distributions, even those not apparent from Fig. 10. The 
TE and LE line sources are particularly well defined. 
The images at and beyond ±18 in. are model-sideplate 
noises and/or sideplate reflections. There are apparent 
phantom images, particularly aft of the TE and around 
the edges of the scanning plane. These are addressed 
below. 

Figures 12 and 13 correspond to Figs. 10 and 11, 
except that DR processing is used for beamforming, 
Eq. (9), and for DAMAS, Eq. (21). It is seen that 
although the beamforming contours differ significantly, 
the source distributions essentially match. The 
exception is that the DR results appear to produce 
cleaner DAMAS results, with much of the phantom 
images removed. That is, the apparent source 
distributions over regions away from surfaces where no 
“real” sources are likely to exist are significantly 
diminished. Also removed with DR is an “apparent” 
LE noise source distribution from the result of 3.15 
kHz. Considering that the present STD method results 
are to some degree contaminated with turbulence 
buffeting microphone self noise, the DR results are 
considered more correct. (This conclusion is supported 
by evidence evaluated in Ref. 5.) 

Also present in Figs. 11 and 13 are edge effects as 
are found in and discussed for Figs. 7 and 8. The edge 
effects can be readily eliminated by expanding the 


scanning frame beyond the regions of strong sources, 
thereby reducing the edge amplitudes and thus any 
potential influence on the regions of interest. This has 
been verified but this is not shown here, as the edge 
effect’s presence in Figs. 11 and 13 is instructive. For 
example, an area where the edge effect appears to 
negatively affect DAMAS results is the sideplate 
region on the left side near the LE (chordwise location 
26 in. and spanwise location -21 in.). The strong array 
responses (Figs. 10 and 12) at that location are not 
correspondingly represented by the DAMAS source 
distributions in that region. Instead, DAMAS puts 
strong sources along the scanning plane edge and the 
LE corner to explain the array response. (Note that it is 
well recognized that the array response over such a 
corner location may well be influenced by reflected 
(and thus correlated) noise sources, whereas the 
DAMAS modeling is based on an equivalent 
statistically independent source distribution. The edge 
effect is unrelated to this modeling/reality physical 
difference. Such reflections undoubtedly cause 
strengthening or weakening and/or shifting of apparent 
sources, but it would not cause source concentration 
along the edges.) Still, even with the scanning plane 
edge effect, away from the edges the TE and LE noise 
source regions are unaffected and the following noise 
spectra serve to verify this. 

Previously 4 , TE noise spectra were determined 
from amplitudes of the array response at the center of 
the TE, along with a transfer function based on an 
assumed line source distribution. Also, corresponding 
spectra from the LE noise region were determined to 
show grit-related LE noise, which due to beamwidth 
characteristics were contaminated by TE noise at low 
frequencies. Figure 14 shows one-third octave spectra 
(per foot) curves from Fig. 13 of Ref. 4 for the test 
conditions corresponding to Figs. 10-13 above. These 
spectra are compared to spectra of TE noise and LE 
noise determined from DAMAS using both STD and 
DR methods. These results are determined by simply 
summing the pressure-squared values of each grid 
point within the rectangular box region surrounding the 
TE and LE regions shown superimposed in Figs. 11 
and 13. The region’s spanwise length is 2.5 feet. The 
sums are divided by 2.5 to put the spectral results on a 
per-foot basis. The spectral comparisons are quite 
good and serve as a strong validation for the different 
analyses. Where low-frequency results of DAMAS are 
not plotted, the integration regions lacked contributions 
(not surprising with the very large beamwidths B). 
The spectra are seen to agree well with the results of 
Ref. 4 over parts of the spectra where each source is 
dominant. Of course in the spectra shown from Ref. 4, 
as in the beamformed solutions of Figs. 10 and 12, the 
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Chordwise location (in.) Chordwise location (in.) Chordwise locations (in.) Chordwise locations (in.) 




FIGURE 10 SADA response contours for Shaded STD FIGURE 11 DAMAS results for Shaded STD processing - TE 

processing - TE and LE noise test. anc l |_(= noise test. Corresponding to Fig. 10. 
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FIGURE 12 SADA response contours for Shaded DR FIGURE 13 DAMAS results for Shaded DR processing - TE 

processing - TE and LE noise test. and LE noise test. Corresponding to Fig. 12. 
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Frequency (Hz) 


FIGURE 14 Comparison of one-third octave spectra from TE 
and LE noise measurements from Ref. 4 and reprocessing by 
DAMAS. 


TE and LE noise region amplitudes are contaminated 
by the mutual influence of the other source over 
different parts of the spectra. The present DAMAS 
results exclude such interference. 

Slat Noise Test. The slat configuration 3 tested in 
the QFF is achieved by removing the flap, filling the 
flap cove (as for the TE noise test above), removing the 
grit boundary layer trip at the LE, tilting the airfoil 
main element to 26° from vertical, mounting the slat, 
and setting the slat angle and gap. The large 26° angle 
is required to obtain proper aerodynamics about the slat 
and LE region 3 . Figure 15 shows the airfoil / slat in the 
deflected open jet, with the SADA positioned 
at0=lO7°. 



FIGURE 15 Sketch of test set up for slat noise test 
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FIGURE 16 Slat noise test, SADA Shaded DR processing, 
/ 1 / 3=20 khlz and Ax/fi=0.083. 



Figure 16 shows the array output and 
corresponding DAMAS result for / 1/3 = 20 kHz using 
shaded DR processing over the scanning plane through 
the airfoil. The distributed slat noise is seen to be well 
identified. There are higher levels toward the left side 
of the slat, likely due to a model /model mount 
irregularity. The aforementioned scanning plane edge 
effect is seen around the edge of the DAMAS 
presentation, and it likely has a mild impact on the 
source definition details at this left side. Away from the 
edge, the slat noise is generally uniform. The 
amplitude of the slat noise is determined by summing 
across the span within the integration box shown. For 
this one-third octave band, the DAMAS level of 57.9 
dB (per foot) compares with 59.1 dB found by Ref. 3, 
which used an approximate procedure involving the 
array output at the slat center and a derived transfer 
function (similar to the procedure of Ref. 4 for TE 
noise). 
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FIGURE 17 Flap noise test, SADA Shaded DR processing, 
/ 1 / 3=20 khlz and Ax/fi=0.083. 


Flap Edge and Flap Cove Noise Test. The flap 
edge noise test configuration is shown in Fig. 6(a) for 
the SADA. The flat edge flap test condition of 29° flap 
angle and M= 0.11 is reexamined for a one-third 
octave frequency band of / 1/3 = 20 kHz. Figure 17 
shows beamforming contours and DAMAS results for 
shaded DR processing over the scanning plane placed 
through the airfoil chordline. The DAMAS results 
appear to successfully isolate the flap edge noise from 
substantial flap cove noise. By using a similar 
rectangular integration region, shown in Figs. 7 and 8, 
one finds a level of 44.6 dB for the flap edge noise. 
This compares to 47.5 dB for a spectrum level 
determined for this flap edge noise case in Ref. 2. In 
that spectrum, this frequency corresponds to a localized 
spectral hump. These present results show that the 
cove noise contributed to the higher level measured. 
DAMAS is seen to allow one to readily separate and 
quantify these cove and flap edge contributions, where 
previously this was not possible. 


CONCLUSION 

DAMAS represents a radical step in array 
processing capabilities. It can replace traditional 
presentations of array results and make the array a 
much more powerful measurement tool than is 
presently the case. The DAMAS equation AX = Y is a 
unique equation that relates a classical beamformed 
array result Y with the source distribution X. The 
sources are taken as distributions of statistically 
independent noise radiators, as does traditional array 
processing/integration analysis. DAMAS does not add 
any additional assumption to the analysis. It merely 
extracts the array characteristics from the source 
definition presentation. The iterative solution for X is 
found to be robust and accurate. Numerical application 
examples show that the actual rate and accuracy at 
which solutions converge depend on chosen spatial 
resolution and evaluation region sizes compared to the 
array beamwidth. Experimental archival data from a 
variety of prior studies are used to validate DAMAS 
quantitatively. The same algorithm is found to be 
equally adept with flap edge/cove, trailing edge, 
leading edge, slat, and calibration noise sources. 
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